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ABSTRACT 



o 



Aims. The role of galaxy mergers in massive galaxy evolution, and in particular to mass assembly and size growth, remains an open question. In 

r ^ this paper we measure the merger fraction and rate, both minor and major, of massive early-type galaxies (M^ > 10^^ M©) in the COSMOS field, 
and study their role in mass and size evolution. 
^ Methods. We use the 30-band photometric catalogue in COSMOS, complemented with the spectroscopy of the zCOSMOS survey, to define close 
^ I pairs with a separation on the sky plane \0h~^ kpc < r^ < 30h~^ kpc and a relative velocity Av < 500 km s~^ in redshift space. We measure both 
major (stellar mass ratio = M^ 2/^*,i > 1/4) and minor (1/10 < < 1/4) merger fractions of massive galaxies, and study their dependence on 
redshift and on morphology (early types vs spirals). 

Results. The merger fraction and rate of massive galaxies evolves as a power-law (1 + z)", with major mergers increasing with redshift, nuu = 1-4, 
^ and minor mergers showing little evolution, n^nm ~ 0. When split by their morphology, the minor merger fraction for early types is higher by a 
, ^ , factor of three than that for spirals, and both are nearly constant with redshift. The fraction of major mergers for massive spirals evolve faster 

(^MM^ ~ 4) than for early types (n^^ = 1.8). 
^-H Conclusions. Our results show that massive early-type galaxies have undergone 0.89 mergers (0.43 major and 0.46 minor) since z ~ 1, leading 
^ ' to a mass growth of ~ 30%. We find that > 1/10 mergers can explain ~ 55% of the observed size evolution of these galaxies since z ~ 1. 
■^j- ' Another ~ 20% is due to the progenitor bias (younger galaxies are more extended) and we estimate that very minor mergers (pi < 1/10) could 
contribute with an extra ~ 20%. The remaining ~ 5% should come from other processes (e.g., adiabatic expansion or observational eff'ects). This 
\^ picture also reproduces the mass growth and velocity dispersion evolution of these galaxies. We conclude from these results, and after exploring 
all the possible uncertainties in our picture, that merging is the main contributor to the size evolution of massive ETGs at z < 1, accounting for 
~ 50 - 75% of that evolution in the last 8 Gyr. Nearly half of the evolution due to mergers is related to minor (pi < 1/4) events. 

Key words. Galaxies: elliptical and lenticular, cD — Galaxies: evolution — Galaxies: interactions 

> 1. Introduction 120061: iPerez-Gonzalez et aP l2008t IPozzetti et alJ lIoTol) . These 
? properties result from several physical mechanisms for which 

p><; _ The history of mass assembly is a major component of the jj necessary to evaluate the relative impact. In this paper we 

H . galaxy formation and evolution scenario. The evolution in the examine the contribution of major and minor mergers to the 

. >^ . number of galaxies of a given mass, as well as the size and j^^ss growth and size evolution of massive early-type galaxies 

shapes of galaxies building the Hubble sequence, provides (eTGs), based on new measurements of the pair fraction from 

strong input to this scenario. The optical colour - magnitude di- (he COSMO^] (Cosmo logical Evo lution Survey, IScoville et al.1 

agram of local galaxies shows two distinct populations: the red ,2007) and zCOSMOSEI dLillv et alJl2007h surveys, 

sequence, consisting primarily of old, spheroid-dominated, qui- ^he number density of red massive g alaxies with ^ 

escent galaxies, and the blue cloud, formed primarily by spi- jqU is roughly constant since z ~ 0.8 (iPozzetti et alJl20T(l 

ral and irregular star-forming galaxies (e.g., Strat eva et al. || 2001| ; ^nd references therein), with major mergers (mass or luminosity 

| Baldryetal. || 2004 ). This bimodahty has been traced at increas- ^atio higher than 1/4) c ommon enough to explain their number 

ingly higher redshifts (e.g., | llbert et al. || 201Q) , showing that the evolution since z = 1 dEliche-Moral et alJl201(i iRobaina et al.1 

most massive galaxies were the first to populate the red sequence [2010; Qe sch et alJl2010h . However, anddespite that they seem 

as a result of the so-called "downsizing" (e.g., | Bundy et al. | -^ead" since z ~ 0.8, two observational facts rule out the pas- 



* Based on observations made at the European Southern Observatory ' http://cosmos.astro.caltech.edu/ 
(ESO) Very Large Telescope (VLT) under Large Program 175.A-0839. ^ http://www.astro.phys.ethz.ch/zCOSMOS/ 
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sive evolution of these massive galaxies after they have reached 
the red sequence: the presence of Recent Star Formation (RSF) 
episodes and their size evolution. In the former, the study of red 
sequence galaxies in the A/^f/V-optical colour vs magnitude dia- 
gram reveals that -30% have undergone RS F, as seen fro m their 
blue NUV - r colours, both locall y (Kavir aj et al .1120071) and at 
higher redshifts (z ~ 0.6, Kaviraj et al. 2 0111). This RSF typically 
involves 5 - 15% of the galaxy stellar mass (IScarlata e t al. 2007; 
iKaviraj et al.l l2008l 1201 ll) . Some authors suggest that minor 
mergers, i.e., the merger of a massive red sequence galaxy with 
a less massive (mass or luminosity ratio lower than 1/4), gas- 
rich satellite, could explain the observed properties of galaxies 
with RSF (Kaviraj et al. 2009 ; Fernandez:j3ntiver os et al. 2011; 
lDesaietal.ll201lh . 

Regarding size evolution, it is now well established that mas- 
sive ETGs have, on average, lower effective radius (re) at high 
redshift than locally, bein g_~ 2 and ~ 4 times smaller at z ~ 1 



and z ~ 2, respectively (Daddi et al.*'2005*: 'Trui illq et alJ 
| 2007L[20 11: Buitrago et al. 2008; van Dokkum et al.l2008L 
'van derWel et al. 2008; Toftetal. 2009; Williams et al. 



2006 



201C 



201C 



Newman et al. 2010, 2012; Damjanov et al. 201 1; Weinzirl et al. 
20111: ICassata et al.1 2011. but see ISaracco et al.1 2010 and 
Valentinuzzi et al. n2010b for a different point of view). ETGs 



as compa ct as observed at h i gh redshifts are rare in the local 
unive rse fTruiillo et al.1 120091: Q^lor et al.l [2Q1Q|: ICassata eTaP 
I2OIII) . suggesting that they must evolve since z ~ 2 to the 
present. It has been proposed that high redshift compact galaxies 
are the cores of present day ellipticals, and that they increased 
their size by adding stellar mass in the outskirts of the galaxy 
(iBezanson et al.|[2009l: iHopkins et al. l|2P09a'; 'van Dokkum et al. 
I201(lh . Several studies suggest that merging, especi ally the mi- 
nor one, could explain the observed size evolution (iNaab et al 



20091: IBezanson et al.||2009|; iHopkins et an2010bl: IShankar et al 
lOseret al. 20 if^ 



pser et al. 201^, while other processes, as adiabatic ex- 
pansion due to AGNs or to the passive evolution of the stellar 
population, should have a rnild ro l e at z < 1 (Fan et al.l l2010l: 
iRagone-Figueroa & Granatoll201 ll: iTruiillo et al.ll201ll) . In ad- 
dition, a significan t fractio n of local ellipticals present signs of 
recent interactions dvanDokkum 2005; Tal et al. 2009). 

While minor mergers are expected to contribute significantly 
to the evolution of massive ETGs, there is no direct observational 
measurement of their contribution yet. As a first attempt, the mi- 
nor merger fraction of the global population of Lb > galax- 
ies in VVDS-Dee|£l (VIMQS VLT Deep S pectroscopic Survey Ban 
iLe Feyre et al.ll2005b has been studied by iLopez-Sanjuan et al.l 
LSll hereafter), who show that minor mergers are quite 
common, that th eir importance decrease with redshift (see also 
iLotz et al]|201ll) . and that participate to about 25% of the mass 
grov ^th by merging o f such g alaxies. Focusing on ma s sive ga lax- 
ies. IWilliams et al.1 (l201ll) . iMarmol-Queralto et al.1 (l2012h . or 
^Newm an et al.l (12012 ) study their total (major -\- minor) merger 
fraction to z ~ 2, finding that it is nearly constant with redshift. 
In this paper we present the detailed merger history, both minor 
and major, of massive (M* > 10^^ M©) ETGs since z ~ 1 using 
close pair statistics in the COSMOS field, and use it to infer the 
role of major and minor mergers in the mass assembly and in the 
size evolution of these systems in the last ~ 8 Gyr. 

The paper is organized as follow. In Sect. [2] we present our 
photometric catalogue in the COSMOS field, while in Sect. [3] 
we review the methodology used to measure close pairs merger 
fractions when photometric redshifts are used. We present our 
merger fractions of massive galaxies in Sect. IH and the inferred 



http ://cesam. oamp . fr/vvdsproj ect/vvds . htm 



merger rates for ETGs in Sect. [5] The role of mergers in the 
mass assembly and in the size evolution of massive galaxies 
is discussed in Sect. [6l and in Sect. |7] we present our conclu- 
sions. Throughout this paper we use a standard cosmology with 
= 0.3, Qa = 0.7, Ho = 100/z Km s"^ Mpc"^ and h = 0.7. 
Magnitudes are given in the AB system. 

2. The COSMOS photometric catalogue 

We use the COSMOS catalogue with photometric r edshifts de- 
rived from 30 broad and medium bands described in lllbert et al.l 
(2009) and iCapak et al.1 (120071) . version 1.8. We restrict our- 
selves to objects with t < 25. The detection completeness at 
this limit is higher than 90% (ICapak et al.ll2007l) . In order to ob- 
tain accurate colours, all the images were degraded to the same 
point spread function (PSF) of 1.5''. At ~ 25, the rms accu- 
racy of the photometric redshifts (Zphot) at z < 1 is ~ 0.04 in 
(^spec - ^phot)/(l + ^spec), where Zsnec IS the spectroscopic redshift 
of the sources (Fig. 9 in lllbert et al.ll2009 ). At z > 1 the quality 
of the photometric redshifts quickly deteriorates. Additionally, 
and because we are interested on minor companions, we require 
a detection in the band to ensure that the stellar mass esti- 
mates are reliable, thus we add the constraint < 24. 

Stellar masses of the photometric catal ogue have been de - 
rived following the same approach than in lllbert et al.l (2010). 
We used stellar populatio n synthesis models to convert liiminos - 
ity into stellar mass (^e.g.. lBell et"aDl2003l: IFontana et al.ll2004l) . 
The stellar mass is the factor needed to rescale the best-fit tem- 
plate (normalized at one solar mass) for the intrinsic luminosi- 
ties. The Spectral Energy Distribution (SED) templates were 
generate d with the stellar populati on synthesis package devel- 
oped by iBruzual & CharlotI (I2003L BC Q3). We assumed a uni- 
versal initial mass function (IMF) from' Chabried (12003 ^ and an 
exponentially declining star formation rate, S FR oc e^^^ (r in 
the range 0.1 Gyr to 30 Gyr). The SEDs were generated for a 
grid of 51 ages (in the range 0.1 Gyr to 14 5 Gyr). Dust extinc- 
tion was applied to the templates using the lCalzetti et al.l (l2000l) 
law, with E(B - V) in the range to 0.5. We used models with 
two diff^erent metalli cities. Following l^ntana et al. ( 2006) and 
IPozzetti et al.l (l2007l) . we imposed the prior E(B - V) < 0.15 if 
age/r > 4 (a significant extinction is only allowed for galaxies 
with a high S FR). The stellar masses derive d in this way have 
a systematic uncertainty of -0.3 dex (^e.g.. IPozzetti et al.ll2007l: 
Barro et al. 2011). 

We supplement the previous photometric catalogue with the 
spectroscopic information from zCOSMOS survey, a large spec- 
troscopic redshift survey in the central area of the COSMOS 
field. In this analysis we use the final release of the bright part 
of this survey, called the zCOSMOS-bright 20k sample. This is 
a pure magnitude selected sample with Iab < 22.5. For a de- 
tailed description and releva nt results of the pr e vious 10k re- 
lease, see Lillv et al. ( 2009); iTasca et al.l (l2009l) : iPozzetti et aP 
( 2010) or Pen g et al.l (12010). A total of 20604 gala xies have been 
observed with the VIMOS spectrograph ( Le Fevre et al.ll20()3h 
in multi-slit mode, and the data have been processed using the 
VIPGI data processing pipeline ( Scodeggio et al.|[2QQ5b . A spec- 
troscopic flag has been assigned to each galaxy provi ding an e s- 
timate of the robustness of the redshift measurement (iLilly et afl 
2007). If a redshift has been measured, the corresponding spec- 
troscopic flag value can be 1, 2, 3, 4 or 9. Flag = 1 meaning 
that the redshift is 70% secure and flag = 4 that the redshift is 
~ 99% secure. Flag = 9 means that the redshift measurement 
relies on one single narrow emission line (O ll or Ha mainly). 
The information about the consistency between photometric and 
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Fig. 1. Stellar mass as a function of redshift in the COSMOS 
field. Red dots are principal galaxies (M^ > 10^^ M©) with 
Zphot in the zCOSMOS area, blue dots are companion galaxies 
(M^ > 10^^ Mo) with Zphot in the COSMOS area, and black dots 
are the red galaxies (NUV-r^ > 3.5) with Zphot in the COSMOS 
area. We only show a random 15% of the total populations for 
visualisation purposes. Green squares mark those galaxies in 
previous populations with a spectroscopic resdhift. The vertical 
lines marks the lower and upper redshift in our study, while the 
horizontal ones the mass selection of the principal (solid) and 
the companion (dashed) samples. 



spectroscopic redshifts has also been included as a decimal in 
the spectroscopic flag. In this study we select the highest reli- 
able redshifts, i.e. with confidence class 4.5, 4.4, 3.5, 3.4, 9.5, 
9.3, and 2.5. This flag selection ensures that 99% of redshifts 
are believed to be reliable based on duplicate objects (iLilly et al.l 
[2009). 

Our final COSMOS catalogue comprises 134028 galaxies at 
0.1 < z < 1.1, our range of interest (see Sect. 12. lb . Nearly 35% 
of the galaxies with < 22.5 have a high reliable spectroscopic 
redshift. For consistency and to avoid systematics, we always 
use the stellar masses and other derived quantities from the pho- 
tometric catalogue. We checked that the dispersion when com- 
paring stellar masses from Zphot and Zspec is ~ 0.15 dex, lower 
than the typical error in the measured stella r masses (- 0.3 dex) 
Thanks to the methodology developed in iLopez-Sanjuan et aD 
( 2010a) we are able to obtain reliable merger fractions from pho- 
tometric catalogues under some quality conditions (Sect. [3]). We 
check that the COSMOS catalogue is adequate for our purposes 
in Sects. [I2]and[33l 

2.1. Definition of ttie mass-selected samples 

We define two samples selected in stellar mass. The first one 
comprises 2047 principal massive galaxies in the zCOSMOS 
area, where spectrosc opic inforniation is available, with > 
10^1 Mo (M^ > M! . lllbertetal.ll2010l) at 0.1 < z < 1.1. The 
second sample comprises the 23992 companion galaxies with 
> 10^^ Mo in the full COSMOS area and in the same red- 
shift range. The mass limit of the compani on sample ensures 
completeness for red galaxies up to z ~ 0.9 (iDrory et al.ll2009l: 
lllbert et al.ll2QTQh . Because of that, we set Zup = 0.9 as the upper 
redshift in our study, while Zdown = 0.2 to probe enough cosmo- 
logical volume. However, our methodology takes into account 
the photometric redshift errors (see Sect. [51 for details), so we 
must include in the samples not only the sources with z < Zup, 



Fig. 2. as a function of redshift in the mass-selected sample, 
from M^ > 10^1 Mo (thiner line) to 10^^ Mo < M^ < 10^^-^ Mo 
(thicker line) galaxies in bins of 0.2 dex. The black solid line 
marks the photometric errors of blue galaxies in the lower mass 
bin, while the dashed line is for red galaxies in the same mass 
bin. The vertical line marks the higher redshift in our samples, 
^max = 1.1. The horizontal line marks the median for low- 
mass galaxies atl<z<l.l,A^ = 0.015. 



but also those sources with z - 2(Tphot < ^up in order to ensure 
completeness in redshift space. Because of this, we set the max- 
imum and minimum redshift of the catalogues to Zmin = 0.1 and 
Zmax = 1.1. We show the mass distribution of our samples as a 
function of z in Fig.[T] and we assume our samples as as volume- 
limited mass-selected in the following. 

Our final goal is to measure the merger fraction and rate of 
massive ETGs, but the our principal sample comprises ETGs, 
spirals and irregulars. We segregate morphologically our princi- 
pal sample thanks to the morphological classification defined in 
iTasca et al.l (120091) . Their method use as morphological indica- 
tor the distance of the galaxies in the multi- space C - A - G 
(Concentration, Asymmetry and Gini coefficient) to the posi- 
tion in this space of a training sample of -500 eye-ball clas- 
sified galaxies. These morphological indices were measured in 
the HST/ACS ima ges of the COSMO S field, taken through the 
w ide F814W filter (iKoekemoer et al.h2007) . We refer the reader 
to lTasca et al.l (l2009l) for further details. The morphological clas- 
sification in the COSMOS field is reliable for galaxies brighter 
than < 24, and all our principal galaxies are brighter than < 
23.5 up to z = 1 . According to the classification presented in 
iTasca et aP (l2009l) our principal sample comprises 1285 (63%) 
ETGs (E/SO) and 632 (31%) spiral galaxies. The remaining 6% 
sources are half irregulars (65 sources) and half massive galax- 
ies without morphological classification (65 sources). The mean 
mass of both ETGs and spirals is similar, M* ~ 10^^-^ Mo. 

2.2. Dependence of photometric errors on stellar mass 

The quality of the photometric redshifts in COS MOS decreases 
for faint objects in the band (lllbert et al .1120091) . In this section 
we study in details how redshift errors depend on the mass of 
the sources, since this imposes limits on our ability of measure 
reliable merger fractions in photometric catalogues (Sect. 13.21) . 
As shown by Jlbert et al.l (I2010l) . we can estimate the photo- 
metric redshift error (cr^ph^t) from the Probability Distribution 
Function of the photometric redshift fit. In Fig. [2] we show 
the median A^^phot = crzphot/(l + ^phot) of galaxies with differ- 
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ent stellar masses, from > 10^^ M© (massive galaxies) to 
10^^ Mq < < 10^^-^ Mo (low-mass galaxies) in bins of 0.2 
dex. 

Massive galaxies are bright in the whole redshift range un- 
der study, thus their photometric errors are small up to z ~ 1, 
A^,phot ~ 0.005. On the other hand, low-mass galaxies are fainter 
at high redshift than their local counterparts, so their their pho- 
tometric errors increase with z and reach A^^phot ~ 0.015 at 
z ~ 1 . We study separately the photometric errors of low-mass 
red and blue galaxies. We take as red galaxies those with SED 
colour NUV - > 3.5, w hile as blue those with NUV - < 
3.5 (see lllbert et al.ll2QTQ[ for details). Blue galaxies also have 
A^,phot ~ 0.015 up to z ~ 1, while red galaxies have higher pho- 
tometric redshift errors, with A^^phot ~ 0.020 at z = 0.95 and 
A^,phot ~ 0.040 at z = 1.05. This different behaviour can be ex- 
plained by the different mass-to-light ratio (M^ /L) of both pop- 
ulations. Faint (/^ ~ 25) blue galaxies, which photometric errors 
are higher, reach masses as low as ~ 10^-^ M© at z ~ 1. On 
the other hand, we are in the detection limit for red galaxies at 
these redshift (these red galaxies have ~ 25, Sect. 12.11) , ex- 
plaining their high photometric redshift errors. In Sect. 13.21 we 
prove that our methodology is able to recover reliable merger 
fractions in COSMOS samples with A^^phot ^ 0.040, as those in 
our study. 

3. Close pairs using photometric redshifts 

The linear distance between two sources can be obtained from 
their projected separation, = 6dA(Zi), and their rest- frame rela- 
tive velocity along the line of sight, Av = c \zj - z/|/(l +Z/), where 
Zi and Zj are the redshift of the principal (more luminous/massive 
galaxy in the pair) and companion galaxy, respectively; is the 
angular separation, in arcsec, of the two galaxies on the sky 
plane; and dA(z) is the angular scale, in kpc/arcsec, at redshift 



With the previous definitions the merger fraction is 



z. Two galaxies are defined as a close pair if r^^^ < 
and Av < Av"^^^. The lower limit in is imposed to avoid 



seeing effects. We used rj 



mm 
P 



lOh-' kpc. 



^P 



30h-^ kpc. 



and Av"^^^ = 500 km s'^ With these const raints 50%-7Q% of 
the selected close pair s will finally merge (iPatt on et al. ' [200QI: 
iPatton & Atfieldl l2Q08l: iLin et al.l l2004l: iBell et al. 2006). The 
PSFo f the COSMOS ground-based images is 1.5'' ( Capak et al." 
l2007h . which corresponds to ~ kpc in our cosmology at 
z ~ 0.9. To ensure well deblended sources and minimize colour 
contamination, we fixed to 10h~^ kpc (6 > 2''). On the other 
hand, we set r^^^ to 30h~^ kpc to ensure reliable merger fractions 
in our study (see Sect. l3.2l for details). 

To compute close pairs we defined a principal and a compan- 
ion sample (Sect. 12. it . The principal sample contains the more 
massive galaxy of the pair, and we looked for those galaxies in 
the companion sample that fulfil the close pair criterion for each 
galaxy of the principal sample. If one principal galaxy has more 
than one close companion, we take each possible pair separately 
(i.e., if the companion galaxies B and C are close to the prin- 
cipal galaxy A, we study the pairs A-B and A-C as indepen- 
dent). In addition, we impose a mass difference between the pair 
members. We denote the ratio between the mass of the principal 
galaxy, M^j, and the companion galaxy, M*,2, as 



^ = 



★,2 



(1) 



and looked for those systems with M^,2 > y^^*,!- We define as 
major companions those close pairs with ju > 1/4, while minor 
companions those with 1/10 < yu < 1/4. 



/m (>//) = 



A^pM 



(2) 



where Ni is the number of sources in the principal sample, and 
A^p the number of principal galaxies with a companion that fulfil 
the close pair criterion for a given yu. This definition applies to 
spectroscopic volume-limited samples. Our samples are volume- 
limited, but com bine spectroscopic and photom etric redshifts. In 
a previous work, iLopez-Sanjuan et aP (120 10a ) developed a sta- 
tistical method to obtain reliable merger fractions from photo- 
metric catalogues. We remind the main points of this method- 
ology below, while we study its limits when applied to our 
COSMOS photometric catalogue in Sect.[T2l 

We use the following procedur e to define a close pair system 
in our photometric catalogue (see iLopez-Sanjuan et aTl l2010aL 
for details): first we search for close spatial companions of a 
principal galaxy, with redshift zi and uncertainty cr^^ , assuming 
that the galaxy is located at zi - 2(T^^ . This defines the maximum 
possible for a given r^^^ in the first instance. If we find a com- 
panion galaxy with redshift Z2 and uncertainty cr^^ in the range 
Tp < r^^^ and with a given mass with respect to the principal 
galaxy, then we study both galaxies in redshift space. For con- 
venience, we assume below that every principal galaxy has, at 
most, one close companion. In this case, our two galaxies could 
be a close pair in the redshift range 



[Z = [Zi - 2cr,^,Zi + 2cr,J [Z2 - 2o-,^,Z2 + 2cr,J. 



(3) 



Because of variation in the range of the function dA(z), 

a sky pair at zi - 2cr^^ might not be a pair at zi + 2cr^^ . We thus 
impose the condition r^^^ < ^p < r^^^ at all z e [z~,z^], and 
redefine this redshift interval if the sky pair condition is not sat- 
isfied at every redshift. After this, our two galaxies define the 
close pair system k in the redshift interval where the 

index k covers all the close pair systems in the sample. 

The next step is to define the number of pairs associated at 
each close pair system k. For this, we suppose in the following 
that a galaxy / in whatever sample is described in redshift space 
by a probability distribution Pi (zi \ T]i), where Zi is the source's 
redshift and rji are the parameters that define the distribution. If 
the source / has a photometric redshift, we assume that 



Pi (Zi I r]i) = Pg (Zi I Zphot,/, cTzphot,/) 



1 I (zi - ZphouiY 

exp 



^f2no■, 



^phot,/ 



while if the source has a spectroscopic redshift 

Pi {Zi I rii) = Pd {Zi I Zspec,/) = ^{Zi - Zspec,/), 



(4) 



(5) 



where d{x) is delta's Dirac function. With this distribution we are 
able to statistically treat all the available information in z space 
and define the number of pairs at redshift zi in system k as 



Vk{Zl) = CkPl{Zl\m) P2(Z2\m)dZ2 



where zi e [z^ , z^], the integration Umits are 
z" =Zi(l-Av™»/c)-Av'"^Vc, 



(6) 



(7) 
(8) 
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the subindex 1 [2] refers to the principal [companion] galaxy in 
k system, and the constant Ck normalizes the function to the total 
number of pairs in the interest range 



2Nl 



X 



Pl(Zl |770dzi + /'2fel^2)dZ2. 



(9) 



Note that Vy^ = if zi < or z\ > z'^. The function Vk (Eq. |[6l) 
tells us how the number of pairs in the system k, N^, are dis- 
tributed in redshift space. The integral in Eq. Q spans those 
redshifts in which the companion galaxy has Av < Av"^^^ for a 
given redshift of the principal galaxy. 

With previous definitions, the number of pairs in the interval 

Zr,l = [ZhZl+l) is 



m,/ 



Pi(Zi\71i)dZi 



(10) 



where the index / spans the redshift intervals defined over the 
redshift range under study. If we integrate over the whole red- 
shift space, Zr = [0, oo], Eq. ([TOb becomes 



(11) 



where Yjk^p is analogous to A^p in Eq. ©. In order to es- 
timate the statistical error of /m./, denoted (Tstat,/, we use the 
jackknife technique (lEfronll 19821) . We compute partial standard 
deviations, Sk, for each system k by taking the diff'erence be- 
tween the measured /m,/ and the same quantity with the ^th 
pair removed for the sample, f^^, such that Sk = fm,i - f^i- 
For a redshift range with A^p systems, the variance is given by 
^stat / = [^^P ~ 1) ^ ?]/A^p- When A^p < 5 we used instead the 
Bayesian approach of ICameronI (|2Q11|) . that provides accurate 
asymmetric confidence intervals in these low statistical cases. 
We check that for A^p > 5 both jackknife and Bayesian methods 
provide similar statistical errors within 10%. 

3.1. Dealing with border effects 

When we search for close companions near to the edges of the 
images it may happen that a fraction of the search volume is 
outside of the surveyed area, lowering artificially the number 
of companions. To deal with this we select as principal galax- 
ies those in the zCOSMOS area, i.e., in the central 1.6 deg^, 
while we select as companions those in the whole photometric 
COSMOS area. This maximize the spectroscopic fraction of the 
principal sample and ensures that we have companions inside all 
the searching volume. 

3.2. Testing the methodology with 20k spectroscopic sources 

Following iLopez-Sanjuan et all (l2QlQah . we test in this sec- 
tion if we are able to obtain reliable merger fractions from our 
COSMOS photometric catalogue. For this, we study the merger 
fraction /m in the zCOSMOS-bright 20k sample. The merger 
fraction in the 10k sa mple is studie d in details by de Rave l et al.l 
( 201 1) and K ampczyk et al.l (1201 ih . We defined /spec as the frac- 
tion of sources on a given sample with spectroscopic redshift. 
The 20k sample has /spec = 1 , while the COSMOS photometric 
catalogue has /pec = 0.34 for < 22.5 galaxies. In this section 
we only use the N = 10542 sources at 0.2 < z < 0.9 with a high 
reliable spectroscopic redshift from the 20k sample. 



To test our method at intermediate /pec, we created synthetic 
catalogues by assigning their measured Zphot and cr^pj^^^ to A^(l - 
/pec) random sources of the 20k sample (we denote this case as 
A = 1 in the following). To explore diff'erent values of A^, we 
assigned to the previous random sources a redshift as drawn for 
a Gaussian distribution with median Zphot and cr^ = (A^ - 1) crj^^^^, 
where A > 1 is the factor by which we increase the initial A^ of 
the sample. In this case, the redshift error of the source is A cr^^^^^ . 
Then, we measure 



syn 



•20k 



(12) 



where is the measured merger fraction in the 20k spectro- 
scopic sample at 0.2 < z < 0.9 without imposing any mass or 
luminosity diff'erence and f^^ is the merger fraction from the 
synthetic samples in the same redshift range. When A > 1, we 
repeated the process ten times and averaged the results. 

We explore several cases with our synthetic catalogues. For 
example, we assume that all sources in the synthetic principal 
catalogue (subindex 1) and in the companion one (subindex 2) 
have a photometric redshift, /pecj = /spec,2 = 0, and that Azj = 
Az,2 = 0.007 (Ai = A2 = 1). We also contemplate more realistic 
cases, as /peci = 0-3 and Azj = 0.007 (Ai = 1) for principals, 
and /spec,2 = and Az,2 = 0.042 (A2 = 6) for companions. We 
found that Sf^ is higher than 10% for r^^^ = 30h~^ kpc close 
pairs for Az,2 ^ 0.05 (A2 > 7) and realistic values of Az,i. We 
checked that \Sf^\ < 10% for Az,2 < 0.04 and r^^^ = 'sO/z"^ 
kpc, justifying the upper limit Az = 0.04 imposed in Sect. 12.21 
For higher r^^^ the method overestimates the merger fraction by 
about 50% in the Az,2 = 0.04 case. Because we are interested on 
faint companions, we set r^^^ = 30h~^ kpc in the following to 
ensure reliable merger fractions. 

On the other hand, we found that the (Tstat of the f^^ is 
~ 5% of the measured value, i.e., two times lower than the es- 
timated I Sfjn I ~ 10%. Because of this, and to ensure reliable 
uncertainties in the merger fractions, we impose a minimum er- 
ror in /m of 10%, and we take as final merger fraction error 
CTf^ = max(0.1/m,crstat). 

In the next section we test further our methodology by 
comparing the merger fraction from a spectroscopic survey 
(/pec = 1) against that in COSMOS from our photometric cat- 
alogue. 



3.3. Comparison with merger fractions in VVDS-Deep: 
cosmic variance effect 

In a previous work in VVDS-Deep, LSll measure the merger 
rate of M| < -20 galaxies with spectroscopic redshifts, where 
M| = Mb + Qz and 2=1.1 accounts for the evolution of the 
luminosity function with redshift, as a function of luminosity 
difference in the 5-band, yUg = Lb,2/Lb,i- As an additional test of 
our methodology, in this section we compare the merger fraction 
in the COSMOS photometric catalogue with that measured by 
LSll down to yUfi = 1/10, reaching the minor merger regime 
in which we are interested on. To minimize the systematic bias, 
we use the same redshift ranges, Zr,i = [0.2,0.65) and Zr,2 = 
[0.65,0.95), close pair definition (r^^^ = 30h~^ kpc), principal 
sample (M| < -20), and companion sample (M^ < -17.5) than 
LSll. We check that the photometric errors are Az < 0.04 up 
to z ~ 0.95 for faint companion galaxies (see Sect. 13.2b . Note 
that LSll use rf"" = 5h-^ kpc, while we take rf"" = lO/z"^ kpc. 
Hence, we recompute the merger fractions in VVDS-Deep for 
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0.18| r 




0.1 0.2 0.3 0.4 0.5 



Fig. 3. Merger fraction of < -20 galaxies as a function of 
luminosity difference in the ^-band, yUg, at z e [0.2, 0.65) (top) 
and z e [0.65,0.95) (bottom) for lOh-^ kpc < < 30h-^ kpc 
close pairs. Diamonds are from present work in COSMOS (pho- 
tometric catalogue) while dots are from VVDS-Deep (LSll, 
spectroscopic catalogue). The black solid lines in both panels 
show the maximum and minimum merger fractions, including 
lo-f^ errors, when we split the COSMOS field in VVDS-Deep 
size subfields (-0.5 deg^). 



^mm _ iQii-i show the merger fractions from COSMOS 

and VVDS-Deep for diff'erent values of j^b in Fig.O 

We find that VVDS-Deep and COSMOS merger fractions 
are in excellent agreement in the first redshift range, while in 
the second redshift range some discrepancies exist, with the 
merger fraction in COSMOS being higher than in VVDS-Deep 
at yu < 1/5. However, both studies are compatible within er- 
ror bars. Note that merger fraction uncertainties in COSMOS 
are ~ 3 times lower than in VVDS-Deep because of the higher 
number of principals in COSMOS. We checked the eff'ect of 
comic variance in this comparison. For that, we split the zCOS- 
MOS area in several VVDS-Deep size (-0.5 deg^) subfields and 
measured the merger fraction in these subfields. The maximum 
and minimum values of /m in these subfields, including IcTf^ er- 
rors, are marked in Fig. [3] with solid lines. We find that, within 
IcTf^, there is a zCOSMOS subfield with merger properties sim- 
ilar to the VVDS-Deep field. Because the zCOSMOS subfields 
are contiguous, this exercise provides a low er limit to the actua l 
cosmic variance m the COSMOS field (e.g., Moster et al.l[20Tl . 
Hence, we conclude that our methodology is able to recover re- 
liable minor merger fractions from photometric samples in the 
COSMOS field. 



4. The merger fraction of massive ETGs in the 
COSMOS field 

The final goal of the present paper is to estimate the role of merg- 
ers (minor and major) in the mass assembly and size evolution 
of massive ETGs. To facilitate future comparison, we present 
first the merger properties of the global massive population in 
Sect. 14.11 Then, we focus in the ETGs population in Sect. 14. 2[ 

The evolution of the merger fraction with redshift up to 
z ~ 1.5 is well parametrized by a power-law function (e.g., 
Le Fevre et al ] |200QI: iLopez-Saniuan et al.ll2009l: |de RaveletalJ 
20S 

/m(z) = /m,0(l+z)^ (13) 

SO we take this parametrization in the following. 

4.1. The merger fraction of the global massive population 

We summarize the minor, major and total merger fractions for 
> 10^1 Mo galaxies in the COSMOS field in Table (Hand 
we show them in Fig. |4l We defined five redshift bins between 
^down = 0.2 and Zup = 0.9 both for minor and major mergers. The 
ranges 0.3 < z < 0.375, 0.7 < z < 0.75 and 825 < z < 85 ar e 
dominated by Large Scale Structures (LSS. lKovac et al.i r2010l), 
so we use these LSS as natural boundaries in our study. This min- 
imizes the impact of LSS in our measurements, since the merger 
fraction depends on environm ent (iLin et al.ll2010l:lde" Ravel et aP 
l2Qlll;lKampczvk et al.""201l). We identify a total of 56.2 major 
mergers and 71.1 minor ones at 0.2 < z < 0.9. Note that the 
number of mergers can be a float number because of the weight- 
ing scheme used in our methodology (Sect.©. We find that 

• The minor merger fraction is nearly constant with redshift, 
/mm ~ 0.051 . The least-squares fit to the minor merger frac- 
tion data is 

= (0.052 ± 0.009)(1 + ^)-o i±o i (14) 

The negative value of the power-law index implies that the 
minor merger fraction decreases slightly with redshift, but it 
is consistent with a null evolution (mmm = 0). This confirms 
the trend found by LS 1 1 for bright galaxies and by lLotz et af] 
( 201 1) for less massive (M* > 10^^ M©) galaxies, and extend 
it to the high mass regime. 

• The major merger rate of massive galaxies increases with 
redshift as 

/mm = (0.019 ± 0.003)(1 + z)'-'^'-\ (15) 

This increase with z contrast with the nearly constant mi- 
nor merger fraction. In Fig. [5] we compare our measurements 
with those from the literature for massive g alaxies and for 
r^^^ ~ 30h-^ kpc close pairs. Ide Ravel et al.l f201 1) use the 
zCOSMOS 10k sample, so their sample is included in ours, 
and r^^^ = 0, r^^^ = 30h~^ kpc spectroscopic close pairs. 
Because they assume a diff'erent inner radius than us, we ap- 
ply a factor 2/3 to their original values (see Sect. [5l for de- 
tails). Both merger fractions are in good agreement, support- 
ing our methodology. Note that o ur uncertaintie s are lo wer 
by a factor of three than those in Ide Ravel et al Jdloni) be- 
cause our princ i pal sa mple is a factor of four larger than 
theirs. IXu et al.l (l2Q12h measure the merger fraction from 
photometric close pairs also in the COSMOS field. They pro- 
vide the fraction of galaxies in close pairs with jj > 1/2.5, 
so we apply a factor 0.7 to obtain the number of close pairs 



6 



C. Lopez-Sanjuan et aL: Mergers and the size evolution of massive ETGs since z ~ 1 





Fig. 4. Major (dots), minor (squares) and total (major + mi- 
nor, triangles) merger fraction of > 10^^ M© galaxies as 
a function of redshift in the COSMOS field. Dashed, solid and 
dott-dashed curves are the least- squares best fit of a power-law 
function, /m oc (I -\- z)"^, to the major (muu = 1.4), minor 
(mmm = -0.1) and total (mm = 0.6) merger fraction data, re- 
spectively. 

(this is the fraction of principal galaxies in their massive sam- 
ple) and a factor 1.6 to estimate the number of yu > 1/4 
systems (the merger fraction depends on yu as /m oc /d\ 
as shown by LSll, and s = -0.95 for mas sive galaxies 
in COSMOS, Sect. O. On the other hand, iBundv et al.1 



(2009) and Bluck et al." (2009) measure the major (// > 1/4) 
merger fraction^ in GOODSQ (Great Observatories Origins 
Deep Survey. iGia valisco et al.l 12004 ) and Palomar/DEEP2 
(IConselice et al.ir2007i) surveys, respectively. These studies 
are also in g ood agreement with our values, with the point at 
z = 0.8 fromlBluck et al.l (I2OO9I) being the only discrepancy. 
The least-squeres fit to all the close pair studies in Fig. \5\ 
yields similar parameters to those from our COSMOS data 
alone, Eq. ([T5l) . 

Regarding morphological studies, iBridge et"aD (1201 Ol ) pro- 
vide the major merger fraction of > 5 x 10^^ M© galax- 
ies in two CF HTLSEI (Canada-Fr ance-Hawaii Telescope 
Legacy Survey. ICoupon et aP l2009l) Deep fields, including 
the COSMOS field. They perform a visual classification of 
the sources, finding 286 merging systems of that mass. We 
cannot compare directly their merger fractio ns with ours 
because of the diff^ere nt methodologies (e.g., iBridge et"aD 
l2010l : lLorzetal.ll201ll) . Thus, we translate they merger rates 
into the expected close pair fraction following the prescrip- 
tions in Sect. \5\ Giving the uncertainties in the merger time 
scales of both methods and the difficulties to assign a pre- 
cise mass ratio yu to the merger ca ndidates in morpholo gical 
studies, the merger fractions from [Bridge et al.l (1201 Ol stars 
in Fig. [5]) are in nice agreement with our results. 
• The fit to the total merger fraction is 



= (0.067 ±0.008)(l-hz) 



0.6±0.3 



(16) 



This evolution is slower than the major merger one, reflect- 
ing the different properties of minor and major mergers. We 
compare our total merger fractions with others in the lit- 
erature in Fig. [6l iMarmol-Queralto et al.l (1201 2h study the 
total merger fraction of massive galaxies by r^^^ = 10h~^ 



Fig. 5. Major (ji > 1/4) merger fraction for > 10^^ M© 
galaxies from r^^^ ~ 30h~^ kpc cl ose pairs. The dots are 
from present work, triangles are form 'de Ra vel et al.l 
the zCOSMOS 10k sample, squares from X u et al.1 (1201 



the COSMOS field, pentagons from Bluck etal. (2009) in the 
Palomar/DEEP2 survey, and diamonds fr om Bundv et al. ( 20091) 
in the GOODS fields. The stars are from IBridge et al.1 (l2010h in 
the CFHTLS by morphological criteria. Some points are slightly 
shifted when needed to avoid overlap. The dashed line is the 
least-squares best fit of a power-law function, /mm (1 -\- z)^'^, 
to the major merger fraction data in the present work. 



The merger fraction depends on the 
-"^•^^ (LSll), so we translate the 



kpc close companions 

search radius as /m c ^ 

merger fractions provided by lMarmol-Queralto et al .l (1 201 
to our search radius. On the other hand. iNewman etaD 
(2012) measure the merger fraction of > 5 x 10^^ M© 
galaxies from r^^^ = 30h~^ kpc close pairs. The values from 
both cl ose pair studies are co nsistent with ours. Also the re- 
sults of I Williams et al.l (12011) suggest a slow/null evolution 
in the total (// > 1/10) merger faction of massive galaxies up 
to z ~ 2. 

Regarding morphological studies, I Jo gee et aP (|2009|) esti- 
mate the total (jj, > 1/10) merger fraction of > 2.5 x 
10^^ Mq galaxies in the GEM^ (Galax y Evolution From 
Morphology And SEDs, iRix et al.l [20041) survey. Their val- 
ues, /m ~ 0.08, are consistent with ours. We also show in the 
merger fraction from Lot z et a l. ( 2011 ) for > 10^^ M© 
galaxies in the AEGISQ (AU-Wavelength Extended Groth 
Strip International Survey. Savis et al.l2007 ) survey. The dif- 
ferent methodologies between these works and ours and the 
different stellar mass regim es probed make dire ct compar- 
isons difficult (see Bridge e t al.ll201Ql: iLotz et al1l2Ql 1 , for a 
review of this topic). In summary, previous work is compat- 
ible with a mild evolution of the total merger fraction, as we 
observe. 



4.2. The merger fraction of ETGs 

We summarize the minor and major merger fractions for both 
ETGs and spiral galaxies in the COSMOS field in Tables[2]and[3l 
respectively, while show them in Fig.|7l We defined five redshift 
bins between Zdown = 0.2 and Zup = 0.9 for ETGs, as for the 
global population, but only three in the case of spirals because 
of the lower number of principal sources. 



^ http://www.stsci.edu/science/goods/ 
^ http://cfht.hawaii.edu/Science/CFHLS/ 



^ http : //www. mpia-hd. mpg . de/GEMS/gems . htm 
^ http://aegis.ucolick.org/ 
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Table 1. Minor, major and total merger fraction of > 10^^ M© galaxies 



Merger fraction 


z = 0.29 


z = 0.46 


z = 0.65 


z = 0.77 


z = 0.86 




0.2 < z < 0.36 


0.36 <z< 0.57 


0.57 <z< 0.73 


0.73 <z< 0.83 


0.83 < z < 0.9 


/mm 


0.031 +0.011 


0.030 + 0.008 


0.035 + 0.007 


0.040 + 0.009 


0.056 + 0.009 


fmm 


0.056 + 0.015 


0.042 + 0.009 


0.047 + 0.008 


0.060 + 0.010 


0.045 + 0.008 


fm 


0.087 + 0.017 


0.072 + 0.011 


0.083 + 0.009 


0.100 + 0.012 


0.101+0.010 



11 > 1/10 



0.20 
0.15 
0.10 
0.05 



0.2 0.4 0.6 0.8 1.0 1.2 

Z 

Fig. 6. Total (major + minor, ji > 1/10) merger fraction as a 
function of redshift. Triangles are from the present work in the 
COSMOS field for ^ 10^ ^ M© galaxies, diamonds are from 
iMdrniol-Queralto et ap (|2012|) for massive galaxies, squares are 
from iNewmanet alf (1201 2h for > 5 x 10^^ M© galaxies, 
crosses are from Josee et al. (I2009h for > 2.5 x 10^^ M© 




;alaxi es by morphological criteria, and dots are from ^ Lotz et al.l 
201 ih for Mi, > 10^^ Mq galaxies by morphological criteria. 
The solid line is the least- squares best fit of a power-law func- 
tion, /m oc (1 -hz)^'^, to the total merger fraction data in the present 
work. 



We assume mmm = in the following for the minor merger 
fraction, as for the global population (Sect. l4.1l) . The mean minor 
merger fraction of ETGs is = 0.060, while /X'""^ = 0.023 for 
spirals. There is therefore a factor of three diff'erence between the 
merger fractions of early type and late type populations. LSll 
also find a similar result when comparing the minor merger frac- 
tion of red and blue bright galaxies. 

On the other hand, the major merger fraction of ETGs is also 
higher than that of spirals by a factor of two. The fit to the major 
merger data yields 



1.8+0.3 



/j^;^ = (0.020 ±0.003)(l+z) 

0.003(1 +zA 



(17) 
(18) 



Because we only have three data points for spirals and of the 
high uncertainty in the first redshift bin, the found value of mmm 
for massive spirals is only tentative. Nevertheless, that the major 
merger fraction of spirals evolves faster than that of ETGs is in 
agreement with previous studies which compare early-types/red 
and late-types/blue galaxies (e.g.,lLin et al Ravel et al.l 

120091; iBundy et aU|2009|;LChou et al.""2011; LSll). 

In summary, the merger fraction of massive (M^ > 10^^ M©) 
ETGs, both major and minor, is higher by a fac tor of 2-3 than 
that of massive spirals (see also iMarmol-Oueralto et al.l l2012l 



0.08 
0.06 
0.04 
0.02 



1 \ \ \ \ \ \ r 



> 10^^ M( 







0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

z 



0.08 



0.06 



1 1 
_ > 10" 


1 1 1 1 1 1 

' , H 


1 1 1 1 1 1 1 1 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Z 

Fig. 7. Major (upper) and minor (lower) merger fractions of 
Mi, > 10^^ Mq galaxies as a function of redshift and morphol- 
ogy. Dots are for ETGs, while squares are for spirals. Dashed 
(solid) lines are the best fit to the early-type (spiral) data, while 
dotted lines are the fits for the global population. 



for a similar result). We estimate the merger rate of ETGs in 
Sect.E 



4.3. Colour properties of companion galaxies 

In this section we attempt to identify the types of galaxies in 
the companion population. As the morphological classification 
is not reliable for all companions because they are faint, we in- 
stead use a colour selection. We split our compani on galaxies 
into red (Wy-r+ > 3.5) and blue (WV-r^ < 3.5. lllbert et aD 
and measure the fraction of red companions (/red) of mas- 
sive galaxies at 0.2 < z < 0.9. 

We find that 62% of the companions of the whole principal 
sample are red, while ~ 38% are blue. Furthermore, the red frac- 
tion remains nearly the same for minor (/red = 60%) and major 
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Table 2. Minor and major merger fraction of ETGs with M^, > 10^^ M© 



Merger fraction 


z = 0.29 


z = 0.48 


z = 0.67 


z = 0.77 


z = 0.86 




0.2 < z < 0.36 


0.36 <z< 0.57 


0.57 <z< 0.73 


0.73 <z< 0.83 


0.83 < z < 0.9 


iw 

J mm 


0.033 + 0.012 
0.056 + 0.016 


0.041 +0.011 
0.054 + 0.011 


0.041 + 0.009 
0.049 + 0.009 


0.053 + 0.012 
0.068 + 0.013 


0.070 + 0.012 
0.070 + 0.011 



Table 3. Minor and major merger fraction of spiral galaxies with > 10^^ Mq 



Merger fraction 


z = 0.39 


z = 0.62 


z = 0.81 




0.2 < z < 0.5 


0.5 < z < 0.7 


0.7 < z < 0.9 


/•spiral 

/mm 

/-spiral 
/mm 


010+^-^22 
025+^-^3^ 


QQ^4+0.011 

0.014 + 0.009 


0.029 + 0.007 
0.029 + 0.009 



(/red = 64%) companions. When we repeat the previous study 
focusing on massive ETGs, we find /red -65%, both for minor 
and major companions. This implies that most of our close pairs 
are "dry" (i.e., red - red). 

5. The merger rate of massive ETGs in the 
COSMOS field 

In this section we estimate the minor (Rmm) and major (Rum) 
merger rate, defined as the number of mergers per galaxy and 
Gyr, of massive ETGs. We remind here the steps to compute the 
merger rate from the merger fraction, focusing first on the major 
merger rate. 

Following lde Ravel et al.l (l2QQ9l) . we define the major merger 
rate as 

^MM = /mm Cp Cm T^^, (19) 

where the factor Cp takes into account the lost companions in 
the inner kpc (iBell et al.l 12006b and the factor Cm is the 

fraction of the observed close pairs that finally merge in a typ- 
ical time scale Tmm- We take Cp = 3/2. The typical merger 
time scale depends on r^^^ and can be estimated by cosmo- 
logical and A^-body simulations. In our case, we compute the 
major merger time s c ale from the cosmological simulations of 
[Kkzbichler & White (^2008), based on the Millennium simula- 
tion (Springel et al. 2005). This major merger time scale refers 
to major mergers (ji > 1/4 in stellar mass), and depends mainly 
on r^^^ and on the stellar mass of the principal galaxy, with 
a weak dependence on redshift in our range of interest (see 
Ide Ravel etal]l2009l for details). Taking log(M^/Mo) = 11.2 
as the average stellar mass of our principal galaxies with a close 
companion, we obtain Tmm = 1.0 ± 0.2 Gyr for r^^^ = 30h~^ 
kpc and Av"^^^ = 500 km s"^. We assume an uncertainty of 
0.2 dex in the average mass of the principal galaxies to esti- 
mate the er ror in Tmm- This time scale already includes the fac- 
tor Cm (see iPatton & Atfieidll2008l: iBundv et al.ll2009l: iLin et al.l 
120 IQl LSll), so we take Cm = 1 in the following. In addit ion, 
LSll show that time scales f rom i Kitzbichler _& White (2008) are 
equivale nt to those from th e A/^-body/hydrodynamical simula- 
tions by iLotz et al.l (l2QlQbh . and that they account properly for 
the observed increase of the merger fraction with r^^^ (see also 



Ide Ravel et al.l l2QQ9l) . We stress that these merger time scales 
have an additional facto r of two uncertainty in t heir normaliza- 
tion (e.g., Hopkins et a 02Q10d:irotz et al.|[2QTTh . 
The minor merger rate is 

^mm — fmm Cp Cm ^mm' (^^) 

where Tmrn = T x Tmm- Following LSll, we take T = 
1 .5 ±0. 1 from the A/^-body/hydrod vnamical si mulations of major 
and minor merg ers performed by iLotz et al.l (20 10b, a, see also 
iLotz et al1l2011l) . As for major mergers, we assume Cp = 3/2 
and Cm = 1 . 

We summarize the merger rates of massive ETGs in Tabled 
and show them in Fig.[8l We parametrize their redshift evolution 
as 

Rm(z)=Rm,0(l+Zr. (21) 

Assuming ^mm = for minor mergers, as for the merger 
fraction (Sect.g!!), we find R^ = 0.060 ± 0.008 Gyr'^ The fit 
of the major merger rate of massive ETGs is 

R^^ = (0.030 ± 0.006) (1 + z)'-^*^-^ Gyr-^ (22) 

Our results imply that the minor merger rate is higher than the 
major merger one at z < 0.5. In addition, the minor and major 
merger rates of massive ETGs are ~ 20% higher than for the 
global population. 

In Fig. [8] we also show the minor and major merger rates of 
red bright galaxies measured by LS 1 1 . We find that red galaxies 
have similar merger rates, both minor and major, as our massive 
ETGs. This suggests that massive red sequence galaxies have 
similar merger properties: nearly 95% of our ETGs are red, while 
the mean mass of the red galaxies in LSI 1 is j-ed ~ 10^^-^ M©, 
a factor of two less massive than our ETGs, ~ 10^^-^ M©. 

The study of the merger properties of red sequence galaxies as a 
function of stellar mass is beyond the scope of this paper and we 
explore this issue in a future work. 

6. The role of mergers in the evolution of massive 
ETGs since z = I 

In this section we use the previous merger rates to estimate 
the number of minor and major mergers per massive (M^ > 
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Table 4. Minor and major merger rate of ETGs with > 10^^ M© 



Merger fraction 


z - 


0.29 


z - 


0.48 


z - 


0.67 


z - 


0.77 


z = 0.86 




0.2 < z < 0.36 


0.36 < 


z < 0.57 


0.57 < 


z < 0.73 


0.73 < 


z < 0.83 


0.83 < z < 0.9 


dET 

mm 


0.049 


+ 0.020 


0.061 


+ 0.021 


0.062 


+ 0.018 


0.080 


+ 0.024 


0.105 + 0.028 


0.056 


+ 0.020 


0.054 


+ 0.016 


0.049 


+ 0.014 


0.068 


+ 0.019 


0.070 + 0.018 



0.14 



I 0.10 h 
^ 0.081- 
^ 0.061- 
^ 0.04 
0.02 




J \ \ \ \ \ \ L 
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^ 0.02 - 
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Z 

Fig. 8. Major (upper) and minor (lower) merger rate of > 
10^^ Mq ETGs as a function of redshift. Filled symbols are from 
the present work, while open ones are from LSll in VVDS- 
Deep for red galaxies. Dashed lines are the best fit to the ETGs 
data, while dotted lines are the fits for the global population. 

10^^ Mq) ETG since z = 1 (Sect. 16. it and the impact of mergers 
in their mass growth (Sect. 16.21) and size evolution (Sect. 16.31) in 
the last ~8 Gyr. 

6.1. Number of minor mergers since z = I 

We can obtain the average number of minor mergers per ETG 
between Z2 and zi < Z2 as 



rz2 t)ET 



dz 



zWoEizY 



(23) 



where E(z) = ^/Q^^^~T~^^^JYl^^ in a flat universe. The defini- 
tion of for major mergers is analogous. Using the merger 
rates in previous section, we obtain A^^^ = 0.89 ± 0.14, with 
A^ET^ = 0.43 ±0.13 and N^^ = 0.46 ± 0.06 between z = 1 and 



z = 0. The number of minor mergers per massive ETGs since 
z = 1 is therefore similar to the number of major ones. Note 
that these values and those reported in the following have an ad- 
ditional factor of two uncertainty due to the uncertainty on the 
merger time scales derived from simulations (Sect. [5]) . 

The number of major mergers per red bright galaxy mea- 
sured by LS 1 1 is = 0.7 + 0.2, higher than our measurement, 
while the number of minor mergers is similar, A/^^ = 0.5 ± 0.2. 
The discrepancy in the major merger case can be explained by 
the evolution of the merger rate in both studies, since LSll as- 
sumed = and we measure = 1.8. 

On the other hand, spiral galaxies have a significantly lower 
number of mergers, A^^^^^^ ~ 0.35, with A^^^^^ ~ 0.15 and 

^mm^ ~ 0.20. We refer the reader to LSll for the discus- 
sion about the role of major and mi nor mergers in the e volu- 
tion of spiral galaxies. In their work, IPozzetti et al.l fcOlQj) find 
that almost all the evolution in the stellar mass function since 
z - 1 is a con sequence of the observed star formation (see also 
IVergani et al.| l2QQ8), and estimate that A/^^ ~ 0-7 mergers since 
z ~ 1 per log(M*/Mo) ~ 10.6 galaxy are needed to explain 
the remaining evolution. Their result is similar to our direct es- 
timation for the global massive population (ETGs -h spirals), 
A^m = 0.75 ±0.14, but they infer A^mm < 0.2. This value is half of 
ours, A^MM = 0.36 ±0.13, pointing out that close pair studies are 
needed to understand accurately the role of major/minor mergers 
in galaxy evolution. 



6.2. Mass assembled through mergers since z = I 

Following LSI 1, we estimate the mass assembled due to mergers 
by weighting the number of mergers in the previous section with 
the average major (//mm) minor merger Gunun) mass ratio, 

SMM) = - 1 = JlMuNHLiO, z) + A^™nA^mm(0, z). (24) 

To obtain the average mass ratios we measure the merger frac- 
tion of massive ETGs at 0.2 < z < 0.9 for diff'erent values of 
yu, from yu = 1/2 to 1/10. Then, we fitted to the data a power- 
law, fm(j^) ^ and used the prescription in LSll to esti- 
mate the average merger mass ratio from the s value. Following 
those steps we find s = -0.95 for massive ETGs in COSMOS, 
while the average merger mass ratios are Jimm - 0.48 and 
y^mm = 0.15, similar to those values reported by LSll. With 
all previous results we obtain that mergers with ji > I / 10 in- 
crease the stellar mass of massive ETGs by 6M^ = 28 ± 8% 
since z = 1. LSI 1 find SM^(1) = 40+ 10% for red bright galax- 
ies in VVDS-Deep, consistent with our measurement within er- 
rors. We note that they use 5-band luminosity as a proxy of 
stellar mass, so their value is an upp er limit due to the lower 
mass-to-light ratio of blue companions. iBluck et al.l(l2QTi]) study 
the major and minor (fi > 1/100) merger fraction of massive 
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galaxies at 1.7 < z < 3 in GNS0 (GOODS NICMOS Survey, 
IConselice et al.l [20TTl) . They extrapolate their results to lower 
redshifts, estimating SM^(l) = 30 ± 25% for yu > 1/10 merg- 
ers. They value is in good agreement with our measurement, but 
its large uncertainty prevents a quantitative comparison. 

The relative contribution of major/minor mergers to our 
inferred mass growth is 75%/25% because the average ma- 
jor merger is three times more massive than the average mi- 
nor on e, as already pointed o ut by LS 1 1 . In their cosmological 
model, iHopkins et al.l (l2010ah predict that the relative contribu- 
tion of major and minor mergers in the spheroids assembly of 
log (Mi, IMq) ~ 1 1.2 galaxies is ~ 80%/20%, in good agreement 
with our observational result. 

On the other hand, several authors have studied luminos- 
ity functions and clustering to constrain the evolution of lu- 
minous red galaxies (LRGs) with redshift, finding that LRGs 
have inc reased their mass SMj, ~ 30%-50% by me rging since 
z = 1 (iBrown et al.1 12007L I2Q08I: ICool et al.1 [2008h . Their re- 
sults are similar to our dir ect estimation, b ut we must take 
this agreement with caution. ^ Tal et al ] (12012) show that LRGs 
have a lack of major companions, excl uding major mergers a s 
an important growth channel (see also [Pe Propris et al.ll2010l) . 
Typically LRGs have L > 3L*, and a low impact of major merg- 
ers in this systems is indeed expected by cosmological models, 
where the contribution of major mergers in galaxy mass assem- 
bly peaks at - MUKho chfar & Silk 2009; Hopkins et al. 2010a; 



ICattaneo et al.ll201 Th /Thus. even if the values of SM^, are sim- 
ilar for LRGs and our massive galaxies, they could have a dif- 
ferent origin. A better approach to estimate indirectly the im- 
pact of mergers in mass growth is to study the evolution of 
massive red galaxies at a fixed number density: because they 
are red (i.e., they have low star formation), their mass is ex- 
pected to growth only by me rging. Followin g this approach, 
Ivan Dokkum et a D (^010) and iBrammer et al ] (1201 Ih estimate 
SMi,(l) ~ 40% for massive galaxi es in th e NEWFIRM Medium- 
Band Surve}0 ( v an Dokkum et al .1120091) . Their result represents 
the integral over all possible // values, so in combination with our 
^M^(l) ~ 30% for jj > 1/10, this would imply that (i) > 1/10 
mergers dominates the mass assembly of massive galaxies since 
z = I and (ii) there is room for an extra SM^, ~ 10% growth due 
to very minor mergers (yu < 1/10). 

6.3. Size growth due to mergers since z = I 

Since the first results of iDaddi et aP (l2005h and iTruiillo et al.l 
("2006), several authors have studied in detail the size evolution 
of massive ETGs with cosmic time. It is now well established 
that ETGs were smaller, on average, than their local counterparts 
of a given stellar mass by a factor of two at z = 1 and of four at 
z = 2 (Sect.[T]). The size evolution is usually parametrized as 



re(0) 



(1+z)-^ 



(25) 



where is the eff'ective radius of the galaxy. Despite of all ob- 
servational eff'orts, the value of a is still in debate, spanning 
the range a = 0.9-1.5 (see references in Sect. [B, as well as 
its dependency on st ellar mass (m assive galaxies e volve faster, 
IWiUiams et al1l2010L or not. .Damianov et al.|[2QTTh . In the fol- 
lowing we assume as fi ducial a value the value reported by 
Ivan der Wei et al.l (l2008l) from a combination of several analy- 



sis, a = 1.2 {Stq - 0.43 at z = 1), with an uncertainty of 0.2 
(dott-dashed line in Fig.O. 

Two main eff'ects could explain the size evolution of ETGs: 
the progenitor bias and genuine size growth. The number 
density of massive (red) galaxies atz = 2 is - 15- 



30% of that in the l ocal universe (e.g.. lArnouts et al.l 120071: 
Pere z-Gonzalez et al.l 120081: IWiUiams et al.l 120101: lllbert et al.l 
2010), and those ETGs that have reached the red sequence 
at later times are systematically more extended than those 
which did it at high redshift. This effect i s called the pro- 
genito r bias and mimic a size gr owth (see Ivan der Wei et al.l 
2009a; Valentinuzzi et al. 2010a,b; Cassata et al. 20111 for fur- 
t her d atails). Both v an de r Wei et al. (2009a) and Sa glia et al.l 
(120101) estimate that the progenitor bias of massive ETGs ac- 
counts for a factor 1.25 (SrQ = 0.8) of the size evolution since 
z = 1, and we assume this value in the following. 

Regarding size growth, several authors have suggested that 
compact galaxies at z ~ 2 are the cores of present day mas- 
sive ellipticals, and that they increase their size by adding stel- 
lar mass in the outski rts of the com pact high redshi ft galaxy 
(Bezanson et al. 2009; Hopkins et al . 2009a; van Dok kum et al.l 
1201 01). The fa ct that the more compact gal axies at z ~ 0.1 
(ITr uiillo et aDi200 9) and z ~ 1 (.Martinez-Manso et al.l [201 ll) 
have similar young ages (~ 1-2 Gyr), combin e d with their 
paucity in the local universe (Truiillo et al.l 120091: iTavlor et al.l 
2010; Cassata etal. 2011) also support the size evolution of 
these systems along cosmic time. Mergers, specially the mi- 
nor ones, have b een proposed to expl ain this evoluti on (e.^ 
iNaab et al.l 120091: iBezanson et al.ll2QQ9l: lHopkins_eTal]l2010b 
Adiabatic expansi on due to AGN activity (F an et al ] |2QTqI) or 
stellar evolution ( Damjanov et al. 2009) could also play a role. 
Thanks to our direct measurements of the minor and major 
merger rate of massive ETGs, we are able to explore the con- 
tribution of mergers to the size growth of these galaxies in the 
last ~ 8 Gyr. 

Theory and simulations show that equal-mass mergers be- 
tween two spheroidal galaxies are less eff'ective in increasing 
the size of ETGs than a major/minor merger with a less dense 
galaxy, both spiral and spheroidal. In the first case the increase 
in size is proportional to the accreted mass, r^ oc M^, with 13=1, 
while in the second case t he index is higher and spans a wide 
range, ~ 1.5 - 2.5 (e.g.. IBezanson et al]|2009l: iHopkins et al.l 
[20 10b). In our case, we estimate J3 for a given ji from the relation 
between the initial (re,i) and the finaleffective radius (re,f ) of an 
ETG in a merger process derived by lFan et aP (l2010l) . 



1 -h 



(26) 



^ http : //www. nottingham. ac . uk/astronomy/gns/ 
^ http : //www. astro . y ale . edu/nmbs/0 verview. html 



^e,i 

wher e e is the slope of the st ellar mass vs size relation. In their 
work. lDamjanov et"aD (|201 ll) find e = 0.5 1 for early-type galax- 
ies in the range 0.2 < z < 2 (see also IWil liams et al. 20101: 
iNewman et al .1120 12b . similar to the 6 = 0.56 from Shen et al] 
(20031) i n SDSS or the e ~ 0.5 e xpected from the Faber- Jackson 
relation (iFaber & Jackson|[l976l) . However, not all the observed 
mergers are between two early-type galaxies. Using colour as 
a proxy for the morphology of our companion galaxies, we 
find that 65% of the mergers are "dry" (red - red), while 35% 
are "mixed" (red - blue), for both major and minor mergers 
(Sect. 14.31) . In the mix ed case we use e = 0.27, a value esti- 
mated from the data of IShen et al.l (l2003l) for late-type galaxies 
in our mass range of interest. Finally, we obtain the JS for a given 
yu as 0.65/5dry + 0.35y gmixed' Using the average mass ratios JImm 
mdJI^ in Sect. 16.21 we find y^MM = 1-30 for major mergers and 
y^mm = 1.65 for minor ones. 



11 



C. Lopez-Sanjuan et aL: Mergers and the size evolution of massive ETGs since z ~ 1 



1.0 
0.9 
0.8 
0.7 h 
^ 0.6 
0.5 h 



/i> 1/10 ' ^ ^ . 

Ii> 1/10 + PB 

All /X + PB 



0.4 1-, , H van der Wel+ 08 
\ \ 



0.0 



0.2 



0.4 0.6 

Z 



0.8 



1.0 



Fig. 9. Effective radius normalized to its local value, dr^, as a 
function of re dshift. The dott-dashed li ne is the observational 
evolution from ' van der Wei et al.1 (l2008h . dr^ = (1 + z)'^'^. The 
solid line is the evolution due to major and minor mergers 
(// > 1/10) expected from our results. The shaded areas in both 
cases mark the 68% confidence interval. The dotted line is the 
expected evolution when the progenitor bias (PB) is taken into 
account. The dashed line is the expected evolution when PB and 
very minor mergers {jd < 1 / 10) are included (see text for details). 

Following Eq. (l24l) , we trace the mass growth of massive 
ETGs with redshift for both minor, ^M^^mmfe), and major merg- 
ers, (5M*,MMb)- Then, we translate these mass growths to a size 
growth with the previous values of 

K fe) = [1 + ^M^,MMfe)]-^^^ X [1 + ^M^,n 



afe)]- 



(27) 



Finally, we estimate the contribution of mergers to the total size 
evolution since z = 1 as 



Are 



1 - drf{\) 
l-^re(l)' 



(28) 



This model yields a size evolution due to mergers of 
drfU) = 0.70 (a = 0.52 + 0.12, solid line in Fig.®. This implies 
that observed major and minor mergers can explain Are -55% 
of the size evolution in massive ETGs since z ~ 1. In the follow- 
ing, all quoted Are have a typical -15% uncertainty due to the 
errors in the merger rates and in the observed size evolution. 

We take into account the progenitor bias by applying a linear 
function to the previous size growth due to mergers (dotted line 
in Fig.O, 



^rf (z) = (l-0.2z)x^r,^(z). 



(29) 



We obtain Sr^^il) = 0.56 (a = 0.84 ± 0.12), thus explaining 
Are ~ 75% of the size evolution with our current observations. 
We note that this value is simi lar to the ^^e (l) = 0.63 estimated 
by the simple model of van der Wei et al.l ([2009 a), which only 
includes the progenitor bias and a merger prescription from cos- 
mological simulations. The remaining Are ~ 25% of the evolu- 
tion should be explained by other physical processes (e.g., very 
minor mergers with ji < 1/10 or adiabatic expansion) or by 
systematic errors in the measurements (e.g., lower merger time 
scales or an overestimation of the size evolution). We explore 
these processes/systematics in the following. 

• Very minor mergers (p < 1/10). Cosmological simu- 
lations find that yu > 1/10 mergers are not the more 
common ones, with the merger history of massive galax- 
ies being dominated by yu < 1/10 mergers (IShankar et al.l 



l2010l:lJimenez et al.l201 lllQser et al.l20T2 !). However, in this 
simulations the mass accretion is dominated by yu > 10 
events due to the low ma ss o f the very minor compan- 
ions. As we show in Sect. 16.21 a mass growth of SM^, ~ 
10% due to very minor mergers since z = 1 is compat- 
ible with the observed n i ass assembly of mas sive galaxies 
(Ivan Dokkum et al ] |2010I: iBrammer et_aD[201ll). This trans- 
lates to A^vm ~ 4 very minor mergers per massive ETG since 
z ~ 1 (we assume that very minor mergers have 1/100 < 
yu < 1/10 and estimate that/Z^j^ ~ 0.025 = 1/40 following 
the prescriptions in Sect. 16.2b . Note that we can increase ar- 
bitrarily the number of very minor mergers by lowering Ji^^, 
but not their contribution to mass growth, which is fixed. We 
checked that the conclusions in this section are independent 

We estimate ySym = 1.85 for very minor mergers, thus ob- 
taining an extra size growth of Are ~ 20% due to mergers, 
Sr^(l) = 0.58 and a = 0.78 ± 0.12 when all yu values are 
taken into account. Hence, mergers since z ~ 1 may explain 
Are ~ 75% of the observed size evolution, while Are -95%, 
with 6rl^(l) = 0.47 and a = 1.1, when the progenitor bias 
is taken into account (dashed line in Fig.O. In this picture, 
nearly half of the evolution due to mergers is related to mi- 
nor (yu < 1/4) events. This result reinforces our conclusion 
that mergers are the main contributors to the size evolution 
of massive ETGs since z = 1 , but observational estimations 
of the very minor merger rate (fi < 1 /10) are needed to con- 
strain t their role. As a first attempt, i Marmol-Queralto et al.l 
(12012 ) find that the merger fraction of massive galaxies at 
z < 1 for yu > 1/100 satellites is two times that of yu > 1/10 
satellites. That suggests A^ym ~ 1, and an additional contribu- 
tion for even smaller satellites (jd < I / 100) could be possible. 
Adiabatic expansion. This will occur in a relaxed system 
that is losing mass. As mass is lost the potential becomes 
shallower, so the system expands into a new stable equilib- 
rium. The amount that a system expands depends on both 
the ejected n iass (M^ ej e ct) and on th e time scale of the pro- 
cess (reject)- iFan et al I (120081 l2010l) suggest adiabatic ex- 
pansion due to quasar activity and/or supernova winds as 
an alternative process to explain the size growth of massive 
early-types, specially at z > 1. These processes occur on 
very short time sca les after the formation of the sphe roid 
(Reject ^ 0-5 Gyr, Ragone-Figueroa & Granatol l201ll) . so 
we expect those galaxies with stellar populations older than 
~ 1 Gyr to be already located in the local stellar mass-size 
relation. This is not supported by observations, in which 
galaxies older than > 3 Gyr at z ~ 1 are still smaller 
than the local ones (see iTrujillo et al.l l201ll for details). 
Interestingly, minor mergers with gas-rich satellites (~ 35% 
of our observed mergers) could trigger recent star formation 
and AGN activity in massive early type s (e.g. jKaviraj et al.l 
l2009l: iFernandez-Qntiveros et al.ll201 ll) . therefore favouring 
some degree of adiabatic expansion and adding an extra size 
growth to the merging process. Devoted N-body simulations 
are needed to explore this topic in details. 
It is also to be noted that the mass loss due to stellar winds 
from the passive evolution of stellar populations in a galaxy 
may lead to adiabatic expansion (Damjanov et al. 2009). 
iRagone-F i gueroa & Gr anato ( 201 1) show that a typical mas- 
sive galaxy is able to eject enough mass due to galactic winds 
to increase its size by a factor of 1.2 in ~ 8 Gyr. This result 
assumes that the potential of the galaxy is not able to retain 
any of the ejected mass, so this could indicate that at most 
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Are ~ 20% of the size evolution since z = I could be ex- 
plained by stellar winds. 

Qverestimation of t he si ze evolution. Results from 
iMartinez-Manso et"an (1201 ih suggest that the photometric 
stellar masses of 'Tru jillo et M (i200 7) are an order of mag- 
nitude higher than those estimated from velocity dispersion 
measurements. This does not erase the size evolution, but 
makes it smaller (massive galaxies are mo re extended than 
less massive ones at a given redshift, e.g.. iDamjanov et al.l 
I201L) . Taking dynamical masses (Mdvm) as a reference in- 
stead of photometric ones. Ivan der Wei et al.l ( 2008) find a = 
0.98 + 0.11, smaller than the a = 1 .20 found by the same au- 
thors from photorn etric studies. The same trend is found by 
ISaglia et al.1 (l2010l) from th e ESO Distant Cluster Surve>(3 
(EDisCSi lWhite et al.l20Q5b galaxies: a ~ 0.65 from dynam- 
ical masses vs a ~ 0.85 from stell ar masses after t he pro gen- 
itor bias is accounted for. Finally. Nevv^man et al.l (1201 Ol) find 
a ~ 0.75 for Mdym > 10^^ M© galaxies. Assuming these 
smaller a values from dynamical masses, major and minor 
mergers account for Are ~ 65% of the size evolution, and all 
the evolution is explained when the progenitor bias and very 
minor mergers are taken into account. 
It is also possible that the extended, low-surface bright- 
ness envelopes of high-z galaxies were missed and their 
were correspondingly underestimated. However, deep obser- 
yations in the near i nfr ared ( optica l res t-frame) from space 
(ISzomom et al.ll2010l: ICassata et al.ll2010|) and from ^round - 
based facilities with adaptive optics (ICarrasco et"al] 120101) 
confirm the compactness of 1 < z < 2.5 massive galaxies. 
On the other hand, also higher values of a than our fiducial 
value a = 1.2 + 0.2 are pr esent in the literature. For ex- 



2008) find a = 1.51 ± 0.04 at z < 2, 



ampl e. iBuitrago et al.l 
while lPamjanov et al.l 

a = 1.5,// > 1/10 mergers would explain Are 



20111) find a = 1.62+0.34. Assuming 

45% of 



the size evolution, while the addition of very minor mergers 
would increase the role of mergers up to Are ~ 65%. In that 
case, the contribution of other processes would increases to 
Are ~ 20%. Thus, even if the size evolution is faster than our 
fiducial a value, mergers would be still the dominant mech- 
anism. 

• Merger time scale. The main uncertainty in our merger 
rates is the assumed merger time scale, which typically 
has a fa ctor of two ii ncertainty in their normalization (e.g., 
[Hopkin s et al.ll2010ch . The Tmm from Kitzbichler & White 
(120081) are ty pically longer than othe rs in the literature 



(^e.g.. iPatton & Atfieldl 120081: II in et al or similar to 

those f rom N-body/hydrodynamical simulations (iLotz et al.l 
I2010bllah . Thus, we expect, if anything, a shorter Tmm, which 
implies a larger role of mergers in size evolution (i.e., higher 
merger rates and numbers of mergers since z ~ 1). In fact, 
a shorter Tmm by a factor of 1.5 is enough to explain the 
observed mass growth and size evolution without the contri- 
bution of very minor mergers. 
• Uncertainties in JS. Equation ([26b . which we used to de- 
rive the values of JS in our model, assumes parabolic orbits 
and dissip ationless (gas-free) nierger s. A bout the first as- 
sumpt ion, iKhochfar & BurkertI (l2006h and IWetzel & Whit3 
^2010) show that most dark matter halos in cosmological 
simulations merge on parabolic orbits. On the other hand, we 
find that ~ 65% of our mergers are dry, but the other -35% 
are mixed and an extra dissipative component is present. In 
these cases simulations suggest that j3 should be higher than 
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Fig. 10. Velocity dispersion normalized to its local value, ^cr^, 
as a function of redshift. The dott-dashed line is the observed 
evolution, Scr^ = (1 + z)^'"^. The shaded area marks the 68% 
confidence interval. The solid line is the evolution due to major 
and minor mergers (yu > 1/10) expected from our results. The 
dotted line is the expected evolution when very minor mergers 
(jd < 1/10) are taken into account. The dashed line is the ex- 
pected evolution when very minor mergers and the progenitor 
bias (PB) are included (see text for details). 



derived from Eq. (l26l) . even reaching J3 - 2.5 (Hop kins et al.l 
l2010b) . This does not change our conclusions because it 
transla tes to a higher si ze evolution due to mergers. In ad- 
dition, 'OseretaLfdloi^) show that the size growth expected 
from Eq. (l26l) is in nice agreement with the growth mea- 
sured in hydrodynamical simulations settled in a cosmologi- 
cal contest. 

In summary, our results suggest that merging is the main con- 
tributor to the size evolution of massive ETGs, accounting for 
Are -50-75% of the observed evolution since z ~ 1. Nearly half 
of the evolution due to mergers is related to minor (ji < 1/4) 
events. 



6.3.1 . Additional constraints from velocity dispersion 
evolution 

In addition to their mass and size, the velocity dispersion (cr*) 
of ETGs evolves with redshift as Scr^ = (I -\- zT. We as- 
su me a = 0.4 ± 0.1 in the following (ScTj, = 1.32 at z = 
1. ICenarro & Truiilloll20Q9l:jCappellari etal.1 120091: ISagfia et al.l 
I2OIOI: Ivan de Sande et al.l |2QTT ). When w e apply our simple 
model using the prescriptions of iFan et al ] (I2OIO) for the evo- 
lution of (T* in merger events, ji > 1/10 mergers are only 
able to explain 15% o f the observed evolution, Scr^(l) = 1.05. 
iHopkins et al ] (I2009bh propose another prescription to trace the 
evolution in cr^ from the evolution in size that takes into account 
the dark matter component of the galaxy. 



y^l/Sr^(z) 



1 -\-y 



(30) 



http://www.mpa-garching.mpg.de/galform/ediscs/index.shtml 



where 7 ~ 1 for M^, ~ 10^^ M© galaxies. Using this prescription, 
the evolution of cr* is faster, but we still explain only ~ 35% of 
the observed evolution, Scr^(l) = 1.10 (solid line in Fig. [TOb . 
The addition of very minor mergers increase the contribution 
to ~ 50%, derail) = 1.16 (dotted line in Fig. [TOI). However, 
a small change of cr^ due to mergers is consistent with the 
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picture from iBernardi et alJ (l2Qlll) . They study in details the 
colour-M^ and colour-cr^ relation of ETGs in SDSS, finding 
that ~ 2 X 10^^ Mq is a transition mass (Mtran) for which the 
curvature of the colour- relation change, while no deviation 
is present in the colour-cr* relation. These authors claim that 
(dry) mergers are the main process in the evolution of those 
ETGs wi th Mj, > Mtran ~ Ml (see also "van der Wei et al." 
2QQ9bl: iLdpez-Saniuan etaP |2Q1Q5 [Oesch et al 2010; 



Eliche-Moral et alJ 



20101: iMendez-Abreu et aP I2012L for 



similar conclusion), as our results also suggest. 

One missing ingredient in the model described in this sec- 
tion is the progenitor bias: new early-types which appeared 
since z ~ 1 are not only more extend ed that previous ones, but 
also have a lower velocity dispersion (Ivan der Wei et al. 2009a). 
Thus, the progenitor bi as also mimic a decrease of cr* with cos- 
mic time. The results of lSaglia et al.l (|2010) suggest that a factor 
of 1 . 1 in the cr^ evolution is due to the progenitor bias. Applying 
this extra evolution as a factor 1 0. Iz to that from mergers (very 
minor ones included), we are able to explain 90% of the increase 
in velocity dispersion, 6cr^^{l) = 1.28 (dashed line in Fig.fTOb. 
Including this, our model is compatible with the observed evo- 
lution and suggests that mergers and the progenitor bias have a 
similar contribution to cr^ evolution, somewhat diff'erent from 
the dominant role of mergers in size evolution. 



6.3.2. Additional constraints from scaling relations 



iNipoti et aP (l2009l) point out that the tightness of the local scal- 
ing laws of ETGs posses an importan t limit to th e growth of 
these systems by (dry) merging (see also lNair et al.|[20 1 1). Using 
these local scaling laws, they conclude that typical present-day 
massive ETGs could not have assembled more than ~ 45% of 
their present stellar mass and grew more than a factor -1.9 in 
size via merging. Even if uncertain, we can extrapolate our ob- 
served trends up to z ~ 2 and compare the inferr ed mass and 
size growths with these upper limits provided by iNipoti et aP 
(^09). We obtain a mass growth by merging (including very 
minor mergers) of SM^, ~ 60% since z = 2, which implies that 
^M^(2)M^(2)/M^(0) ~ 40% of the total mass at z = was 
assembled by merging since z = 2. The size grows by a fac- 
tor of ~ 2 due to this merging in the same cosmic time lapse. 
Therefore, merging seems compatible with the upper limits in 
mass assembly and size growth imposed by the tightness of the 
local scaling laws, although a more complex model is needed to 
fully explore how these laws evolve due to our observed merger 
history. 



6.3.3. Comparison with previous studies 



In a previous work, iTruiillo etal.l(l2011l) use a similar model than 
ours to estimate the number of mergers needed since z ~ 1 to 
explain size evolution if merging is the only process involved. 
They conclude that A^m = 5.0 ± 1.64 mergers with yu = 1/3 
are needed. This number of mergers is higher than our direct 
measurement by a factor of five, A^m = 0-89 ± 0.14 (our average 
merger with ji > 1/10 has /I ~ 1/3). If we take into account 
our estimated very minor mergers, our numbers are ~ 5 and 
/7 ~ 1 / 10. For this value of ji they infer A^^i = 11.20 + 3.66, still 
higher than our estimation. The model of Trujillo et al.l (1201 ll) 
also estimates the mass growth due to mergers since z ~ 1 , which 
is a factor of 3 - 5, also higher than any observational estimation 
or constraint (a factor of ~ 1.4, Sect l6 21) . 



iNewman et al.l (l2012h study the size evolution of red galax- 
ies in the CANDELSM (Cosmic Assembly Near-infrared Deep 
Extragalactic Legacy Survey, Grogin et al. 201 1) survey and the 
role of mergers with jd > 1/10 in this size growth at 0.4 < z < 2. 
Applying a similar model than ours to translate their observed 
total merger fraction to a size growth, they conclude that merg- 
ing can reasonably account for the size evolution observed at 
z < 1 after the progenitor bias is taking into account, while at 
z > 1 mergers are not common enough. Despite they only have 
one merger fraction data point at 0.4 < z < 1 (Fig. O, their 
conclusion is consistent with our more detailed study at z < 1 . 

6.3.4. Expectations from cosmological models 

Several theoretical eff'orts have been conducted to explain the 
size evolution of ETGs. In this section we compare the predicted 
size evolution from cosmological models with our best model, 
which suggests that Are ~ 75% of the evolution in size is due 
to mergers. Are ~ 20% to the progenitor bias and Are ~ 5% to 
other processes (e.g., adiabatic expansion). 

The model of Hopkins et al. (2010b) predicts that, since 
z = 2, un-equal mass mergers explain Are ~ 60% of the ob- 
served size evolution, in agreement with our result. However, 
these authors only track the evolution of compact galaxies since 
z = 2 and do not take into account the possible contribution of 
the progenitor bias, but argue that it should impact their predic- 
tions. In fact, they predict that ~ 45% of the size evolution since 
z = 1 is due to un-equal mass mergers, another ~ 45% is ac- 
counted for by systematics in size measurements and the extra 
~ 10% is due to adiabatic expansion, probably reflecting their 
biased population. 

The model of IShankar et al.l (l201ll) predicts SrQ ~ 0.7 for 
massive galaxies, in agreement with our observational deriva- 
tion due only to mergers (see also iKhochfar & SilP l2006h . 
Interestingly, the evolution increases to ^re ~ 0.5 when individ- 
ual galaxies are tracked along their evolution without any stellar 
mass selection. They predict that ~ 40% of the mass accreted 
by merging in massive galaxies is due to major mergers with 
jd > 1/3. Our best model implies that ~ 47% of mass and size 
growth is due to major mergers with ju > 1/3. The qualitative 
agreement between both works is remarkable. 

On the other hand, Oser et al.l (|2012) find a = lA2(a = 1A4 
for passive galaxies) and a ~ 0.4 by re-simulating with high res- 
olution a set of 40 galaxies with M* > 6.3 x 10^^ M© in a cosmo- 
logical context. They find that the number- averaged merger has 
Ji = 1/16, while the mass-averaged merger has Ji^ = 1/5. From 
our model we estimate Ji^ = 1/3 and Ji = 1/13. We check that 
Ji^ is independent of the assumed number of very minor mergers 
A^vm, while we can vary Ji arbitrarily by changing Nyj^- Thus, 
only the comparison with Ji^ is representative. The predicted 
value is lower than our measurement, but they find a hi gher role 
of major mergers at M^ ~ ipn i Mq (see also IKhochfar & Sil5 



or majc 
[200l;ll 



: ]Hopkins et a02010al:[Cattaneo et al.l201 ll) . with Ji .^ ~ 1/3 
and a big dispersion due to the low statistics (see their Fig. 6). 
Future simulations with higher number of galaxies are needed to 
explore in more details this issue. 

In summary, our result that merging is the main process 
involved in size evolution mostly agrees with simulations, but 
more observational and theoretical studies are needed to under- 
stand the remaining discrepancies. 



http://candels.ucolick.org/About.html 
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7. Conclusions 

We have measured the minor and major merger fraction and rate 
of massive (M^ > 10^^ M©) galaxies from close pairs in the 
COSMOS field, and explored the role of mergers in the mass 
growth and size evolution of massive ETGs since z ~ 1 . 

We find that the merger fraction and rate of massive galaxies 
evolves as a power-law (1 + z)", with no or only small evolution 
of the minor merger rate, ^mm ~ 0, in contrast with the increase 
of the major merger rate, ^mm = 1-4. The total (major + minor) 
merger rate evolves slower than the major one, with rim = 0-6. 
When splitting galaxies according to their HST morphology, the 
minor merger fraction for ETGs is higher by a factor of three 
than that for spirals, and both are nearly constant with redshift. 
The fraction of major mergers for massive spirals evolve faster 
(^MM - 4) than for ETGs (^^t^ = 1.8). 

Our results imply that massive ETGs have undergone 0.89 
mergers (0.43 major and 0.46 minor) since z ~ 1, leading to a 
mass growth of ~ 30% (75%/25% due to major/minor mergers). 
We use a simple model to translate the estimated mass growth 
due to mergers into an effective radius growth. With this model 
we find that ji > 1/10 mergers can explain ~ 55% of the ob- 
served size evolution since z ~ 1. We infer that another ~ 20% is 
due to the progenitor bias (the new ETGs appeared since z = 1 
are more extended than their high-z counterparts) and we es- 
timate that very minor mergers (// < 1/10) could contribute 
with an additional ~ 20%. The remaining -5% could come 
from adiabatic expansion due to stellar winds or from observa- 
tional effects. In addition, our picture also reproduces the mass 
growth and the velocity dispersion evolution of these massive 
ETGs galaxies since z ~ 1 . 

We conclude from these results, and after exploring all the 
possible uncertainties in our model, that merging is the main 
contributor to the size evolution of massive ETGs at z < 1 , ac- 
counting for ~ 50 - 75% of that evolution in the last 8 Gyr. 
Nearly half of the evolution due to mergers is related to minor 
(pi < 1/4) events. 

Studies in larger sky areas are needed to improve the statis- 
tics, especially at lower redshifts when the cosmological vol- 
ume probed is still the main source of uncertainty. We point 
out that a local measurement of the minor merger fraction and 
rate is needed to better constrain its evolution with redshift. 
Understanding the dependency of the minor merger rate on stel- 
lar mass, as well as extending observations to the very minor 
merger regime (fi < 1/10) will be important to further improve 
this picture. In addition, extending the observational work at 
z > 1, when the massive red sequence seems to emerge, will 
be necessary to probe the early epochs of mass assembly. 
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